The Pion Form Factor. 
Where Does It Come From ? 
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Abstract 

We discuss the nonleading, "soft" contribution to the pion form 

' factor at intermediate momentum transfer within operator product 

Q>^ ■ expansion approach. We argue, that the corresponding contribution 

can temporarily simulate the leading twist behavior in the extent 

P^; region of : 3GeV^ < < 35GeV^, where Q^F{Q^) ~ const. 

Such a mechanism, if it is correct, would be an explanation of the 

phenomenological success of the dimensional counting rules (which 

theoretically correspond to the keeping of the asymptotically leading 

ff^ ■ terms only) at available, very modest Q^. The relation with quark 

^ , model calculation is also discussed. 
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1. Introduction 



The investigation of the exclusive asymptotic processes has a long history. 
In early seventies the famous dimensional counting rules were proposed 
The predictions of these rules agree well with experimental data, such as the 
pion and nucleon form factors, large angle elastic scattering cross sections 
and so on. This agreement served as a stimulus for further theoretical in- 
vestigations. The modern approach to exclusive processes was started in the 
late seventies and early eighties |0. We refer to the review papers for 
the detail analysis and discussions. 

However, since 1981, the applicability of the approach 0- |^ at exper- 
imentally accessible momentum transfers was questioned In these 
papers it was demonstrated, that the perturbative, asymptotically leading 
contribution, is much smaller than the nonleading contributions. Similar 
conclusion, supporting this result, came from the different side, from the 
QCD sum rules, where the direct calculation of the form factor has 
been presented at < 3GeV^. This method, unfortunately, by some tech- 
nical reasons, can not directly be applied for the analysis of the form factor 
FniQ"^) (see, however, 0) at larger Q'^ > 3GeV^. Thus, the question " what 
kind of contributions are responsible for the F^^IQ"^) at > 3GeV^ " can 
not be answered within this method. 

At the same time, the information which can be extracted from the dif- 
ferent QCD sum rules, unambiguously shows the asymmetric form of the 
leading twist vr -meson wave function wf, |l^. The application of this wf to 
different amplitudes gives very sensible result at ~ lOGeV"^ and one could 
think that the region 3GeV^ < < lOGeV'^ is the transition zone, where 
the asymptotically leading contribution comes into the game. However, the 
recent papers |jll|,[|l^] do not support this conjecture. Namely, it was shown 



that the inclusion of the intrinsic k'j_- dependence and Sudakov suppression 
leads to the self-consistent calculation, but the obtained magnitude is too 
small (at least there is factor 3 in the description of the vr-meson form factor 
at ~ lOGeV"^) with respect to the data even if (pcz asymmetric function 
is used. 

Now we are ready to formulate the question, which we want to discuss in 
the present letter. 

• If the asymptotically leading contribution can not provide the experi- 
mentally observable absolute values, than how can one explain the very good 
agreement the experimental data with dimensional counting rules [|, which 
supposed to be valid only in the region where the leading terms dominate? 



^ Let us remind, that these rules unambiguously predict the dependence of ampli- 
tudes on dimensional parameters. In particular, Q^Ft^{Q^) ~ constant., Q'^Fp{Q^) ~ 
constant., s'^ ^{iP 7r+n) = f{t/s) .. The experimental data are in a good agreement 
with these predictions in the large region of s,Q'^.&t very modest energy and momentum 
transfer. 
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It is clear, that the possible explanation can not be related to the specific 
amplitude, but instead, it should be connected, somehow, with the nonper- 
turbative wave functions of the light hadrons (vr, p, p...) which enter to the 
formulae for exclusive processes. The analysis of the tt meson form factor, 
presented below supports this idea. 

To anticipate the events we would like to formulate here the result of this 
letter. The very unusual properties of the transverse momentum distribution 
of the nonperturbative vr meson wave function lead to the temporarily 
simulation of the dimensional counting rules by soft mechanism for the 
Ftt{Q'^) at the extent range of intermediate momentum transfer: 3GeV^ < 

< 35GeV^. In this region the soft contribution to Ft,{Q'^)Q'^ does not fall- 
off, as naively one could expect, and we estimate it as F.^{Q'^)Q'^ ^ OAGeV"^. 
The leading twist contribution, after Sudakov suppression, gives, according 
to 0, a little bit less (about 0.2^6^^). 



2. Constraints on the nonperturbative wave function tplk 



±5 ■ 



First of all let me remind some essential definitions and results from [|T^ 



about nonperturbative wf. We define the pion axial wave function in the 
following gauge-invariant way: 

tUq,Mzq,z') = mz)^,j,e''^-^^^'''^ui~z)\niq)) (1) 
= E ^(0|rf(0)7;.75(^^. 5.)"n(0)|7r(g)), 



where Dy=Dy — and = id^j,+gA'^^^ is the covariant derivative. From 
its definition is clear that the set of different vr meson matrix elements defines 
the nonperturbative wave function. The most important (at asymptoticaly 
high Q^) is the part related to longitudinal distribution. The corresponding 
problem was studied in [1T0| , ^. At this paper we are mainly interested in 



transverse distribution part, which is very important at available, not very 
high Q^. We define the mean values of the quark transverse distribution in 
the following way: 

{Q\HlSD,t,fMn{q)) = ^/.g.(-t')"^^|^^(^i")A, (2) 



where transverse vector = (0, t, 0) is perpendicular to the hadron momen- 

2ra-i; 
(2n)!' 



tum = (go, Oj_, qz). The factor '■^.^^.^j" is related to the integration over 



angle in the transverse plane: / (i0(cos0)^"/ J d(j) = {2n — 1)!!/(2?t,)!!. 

We interpret the (fc^) in this equation as a mean value of the quark 
perpendicular momentum. Of course it is different from the naive, gauge 
dependent definition like (0| J7j,759^u|7r(g)), because the physical transverse 
gluon is participant of this definition. However, in QCD the only questions 
which allowed to be asked look like that: "what the magnitude of this gauge 



2 



invariant matrix element?". Our definition (|l|), (|^) does satisfy to tfiis re- 
quirement. We believe that such definition is the useful generalization of the 
transverse momentum conception for the interactive quark system. 

We can not prove that the mean value of the transverse momentum distri- 
bution (A;^), defined above in terms of QCD matrix elements coincides with 
the perpendicular momentum in simple quark model, which is frequently used 
in the phenomenological analysis^. However we would like to make such as- 
sumption in the following section devoted to the phenomenological analysis 
of the pion form factor. At the same time, in this section we prefer to stick 
around QCD with its lanquage of operators, matrix elements and vacuum 
condensates in order to make difference between pure QCD results (without 
any additional assumptions) and phenomenological analysis ( next section), 
where such kind of assumption, motivated by quark model, has been made. 

It turns out, that the calculation of (^i"), defined in (|^), can be reduced 
to the problem of calculating of the mixed vacuum condensates. In particular 

Ckl) = 4 <^'^^"7_^P-^^^ ^ ^ ^ O.lGeV^ ml^O.SGeV^. (3) 



(^1) ..or. (9'G-uG%) 



3K - — ^^-^-5-7, 



where we have introduced the coefficient of nonfactorizability K {K = 1 if 
the factorization would work) of the mixed vacuum condensates of dimension 
seven. This coefficient can be estimated with high enough accuracy and we 
expect K ~ 2.5 ^ 3. 

Few comments are in order. We define the nonperturbative wave func- 
tion ip{k'^,C,), which is Fourier transform of the (j)A{zq, z"^), through its mo- 
ments which can be expressed in terms of nonperturbative vacuum conden- 
sates. As is known, the condensates are defined in such a way, that all gluon's 
and quark's virtualities, smaller than some parameter fi (point of normaliza- 
tion) are hidden in the definition of the "nonperturbative vacuum matrix 
elements" . All virtualities larger than that should be taken into account per- 
turbatively, see |n| for detail discussion of this problem. At the same time. 



there is a relation between moments and condensates, as explained above. 
Thus, all transverse moments are defined in the same way as condensates do. 
This is, actually, the definition of the nonperturbative wave function 
%l){k\,E,)-, through its moments which can be expressed in terms of nonper- 
turbative vacuum condensates. The important consequence of the definition 

In order to prove such identification we need to demonstrate the confinement and 
chiral symmetry breaking , to calculate the meson masses and wave functions and finally, 
to derive an effective low-energy quark model in terms of the fundamental QCD parameters 
(like chiral condensate). We do not have such ambitious purposes at the moment. Instead, 
we assume that such kind of correspondence with intuitive quark model, takes place with 
high enough accuracy. 
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can be formulated in the following way: The nonperturbative wf even for 
sufficiently large /c^ does not behave like l/k\ as is frequently assumed. 

Now we want to formulate some additional constraints on the nonper- 
turbative wave function ip{k\,x) in order to model it. As is known, the 
knowledge of a finite number of moments is not sufficient to completely de- 
termine the iplk^jx); the behavior of the asymptotically distant terms is a 
very important thing as well. We assume, that the 7r-meson fills a finite 
duality interval in the dispersion relation at any n, where n is the order of 
moment: (fcj")- K gives the following constraint: 



1 + e 



For the n = we reproduce the well-known result ||10[ for the wave function: 
(j){C, —>■ ±1) = / dk'j_tp{k'j_, C, — * ±1) ~ (1 — ^^). This constraint is extremely 
important and implies that the kj_ dependence of the ipik"^, ^ = 2x — l) comes 
exclusively in the combination k'j_/{l — at ^ — > ±1. The byproduct 
of this constraint can be formulated as follows. The standard assumption 
on factorizability of the iplk^,^) = ilj{k\)(j){^) does contradict to the very 
general property of the theory formulated above. 

Our last constraint comes from the analysis of the asymptotically distant 
terms. It was argued in that the moments {k"^) at asymptotically large 



n behave in the following way: 

d^Jdki^i^^.oikir^i, 

which means that the =^ ^ij^ " sufficiently large k'j_ 0; the 

S is some input dimensional parameter, which will be expressed in terms of 
{kl). 

We want to model the simplest version of the wave function which meets 
all requirements formulated above. First of all, as was explained, the k'j_ 

dependence comes only in the combination ipijz^) and at sufficiently large 
k'j_ we have to have the — S) -like dependence. Besides that, in order 

to reproduce the noticeable fluctuations of the k1_ (|^), we have to spread out 
the distribution function between fc^ ~ and k"^ ^ S in such a way, that the 
overall area will be the same, but moments should satisfy to this requirement. 
In principle, it can be done in arbitrary way. The simplest way to make the 
wf wider is to put another 6 function at kj_ = 0. With these remarks in 
mind we propose the following "two-hump" nonperturbative wave function 



^From the physical point of view it is clear that such i5 dependence should be somehow 
regularized. It turns out that the large a" corrections are responsible for this regularization 
(p5[, be published ). The problem of regularization is not an essential point at the moment. 
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which meets all requirements discussed above: 

Hki, = mj^, -s)+ B5i^^)Me + 

(^i) ~ 5(^1)2 ^ A = 7/8, A + B = ^ g = l, ^ = y (fc^ ) ~ o.SGeV^^ 

Few comments are in order. First, this is not supposed to be a realistic model. 
But we expect that this function illustrates features which probably 

does exhibit. We put the common factor [gi^"^ — ^) + |] in the front of the 
formula in order to reproduce the light cone 0(^) function with arbitrary 
(^^). For g = it corresponds to the asymptotic wf : = 3/4(1 — ^^). 
For 5( = 1 we reproduce the 0(Ocz = 15/4(1 - ^^)^^ with {^'^) ~ 0.43. 

Before to consider some applications, let me stress, that we discuss in this 
letter some qualitative properties of the wf. Thus, the model wf proposed 
above should be considered as illustrative example only. We try to emphasize 
on the very unusual qualitative feature of the model. We expect, however, 
that this feature will survive in the full picture, after physical regularization 
will be made. 

3. The Pion Form factor. 

The starting point is the famous Drell-Yan formula [|l6l (for modern, 
QCD- motivated employing of this formula, see [Q), where the F^^^q^) is 
expressed in terms of full wave functions: 

ig^3 fc± + X2q±)^BL{xi, (4) 

where = —q]_ is the momentum transfer. In this formula, the bl{xi., k±) 
is the full wave function; the perturbative tail of bl^xi, k±) behaves as 
as/kj_ for large k1_ and should be taken into account explicitly in the calcu- 
lations. This gives the one-gluon-exchange (asymptotically leading) formula 
for the form factor in terms of soft pion wf with removed perturbative tail 

§• 

Let us remind, that the formula (^) takes into account only the valence 
Fock states. Besides that, k\ in this formula is the usual (not covariant) 
perpendicular momentum of the constituents, and not the mean value {kj_) 
defined in QCD, as a gauge invariant object. However we make the assump- 
tion that it is one and the same variable. The physics behind of it can be 
explained in the following way. 

In the formula (^ we effectively take into account some gluons (not all 
of them), which inevitably are participants of our definition of wf. These 
gluons mainly carry the transverse momentum (which anyhow, does not ex- 
ceed QCD scale of order fi ~ IGeV) or small ammount of longitudinal mo- 
mentum. The contributions of the gluons carrying the finite longitudinal 
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momentum fraction are neglected in (^). This is the main assumption. It 
can be justified by the direct calculation ^ of quark-antiquark-gluon (with 
finite momentum fraction) contribution to vr meson form factor at large 
within the standard technique of the operator product expansion. By techni- 
cal reasons the corresponding calculation has not been completed, however it 
was found that the characteristic scale which enters into the game is of order 
IGeV^. Thus, it is very unlikely to expect that these contributions might 
be important at ~ lOGeV'^. The second calculation, which confirms this 
point, comes from the light cone QCD sum rules [T^. This is almost model 
independent calculation demonstrates that the quark-antiquark-gluon (with 
finite momentum fraction) contribution does not exceed 20% at available Q^. 

Thus, we expect, that by taking into account the only "soft" gluon contri- 
bution (hidden in the definition of kj_ (|^)), we catch the main effect. Again, 
there is no proof for that within QCD, and the only argumentation which 
can be delivered now in favor of it, is based on the intuitive picture of quark 
model, where current quark and soft gluons form a constituent quark with 
nonzero mass and with original quantum numbers. No evidence where a 
gluon would play the role of a valent participant with a finite ammount of 
momentum, is found. 

From the viewpoint of the operator product expansion, the assumption 
formulated above, corresponds to summing up a subset of higher-dimension 
power corrections. This subset actually is formed from the infinite number 
of soft gluons and unambiguously singled out by the definition of nonpertur- 
bative wf (||). 

In the following, we preserve the notation "^Biixi, k±) for the nonpertur- 
bative, soft part only. It should not confuse the reader. 

The formula (^ is written in terms of Brodsky and Lepage notations 
the relation to our wave function k±) looks as follows: 

^bl{xu k±) = ^^^^(e, J j^'^BLixu k^) = (5) 

where /tt = 133MeV. There are two, physically different contributions to 
the form factor. We want to separate them and discuss independently. 

The first contribution is the standard one, in a sense that it is determined 
by the k1_- distribution about k"^ = region. It is proportional to the (we 
call the corresponding contribution F^^) and actually depends very strongly 
on regularization procedure. However, at ~ fewGeV^ this contribution 
begins to fall off very rapidly because of the proportionality to c{k\)'^/ with 
some coefficient c, determined by characteristic scale of the wf. Besides that, 
we expect that all Fock states are equally important in this region and thus, 
any predictions are model dependent. 

However, to give some qualitative insight on the behavior of the contribu- 
tion of this kind, we estimate it in the following way. We regularize our S func- 
tions in the simplest way- we shift its argument: 6{^^) ^{j^ — m^), by 
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introducing the new effective phenomenological parameter m^. This param- 
eter effectively describes the kj_- distribution about zero point and somehow 
related to {k'^)- This regularization is still not sufficient for the calculating 
of the integral at g = 0. But we do not need it, because we are not going 
to apply this method for the describing of the form factor in the vicinity of 
q = 0. Nevertheless, we expect, that at ^ only the integral char- 
acteristic does matter. Have these few comments in mind, we estimate the 
parameter from the requirement of coincidence theoretical formula (§) 
with the experimental data at ~ SGeV'^ : F^^c^ ~ OAGeV'^. It gives 
~ (0.226*6^^)^ for our set of parameters @. We display the result of this 
contribution to the form factor g^F^^lq^) in Fig. 1 by small-dashed curve. 
We took into account in this numerical calculation the evolution of the wf 
which determined by parameter glq"^) 

9{q) = 9if^)i^,r^'\ ^?(/i) = l, /i = 0.5GeV, Aqcd = lOOMeK (6) 

Let me stress again: we are not pretending to have made a reliable calcu- 
lation of the form factor here; we displayed this contribution only for the 
illustrative purposes. The main feature of this contribution- it gives very 
reasonable magnitude for the intermediate region about few GeV^ and it 
starts to fall off very quickly at g^. ^ f^'^ {^D- We expect that any rea- 
sonable, well localized wf with the scale ~ {k'j_) leads to the same behavior. 
We will not discuss this standard contribution any more. 

Currently, much more interesting for us, is the interference term, pro- 
portional to AB. It is clear, that with increasing of g^, the contribu- 
tion coming from the overlap of the 6 functions starts to grow. Because 
the 6 functions are well separated from each other on the value of order 
S ~ O-SGeV"^, we expect that this contribution will start to increase at high 
enough g^ ^ 45* ~ AGeV'^. This contribution almost model- independent 
in a sense that it does not depend on the infrared regularization parameter 

which we introduced before, provided that <^ 5*. Let me remind at 
this point, that our parameters 5* and (fc^) are not the new phenomenolog- 
ical input parameters, but have been derived ^ in terms of the vacuum 
characteristics, which are known independently. 

With these remarks in mind, we can explicitly calculate this interference 
contribution: 

(gl) = ^5^(1 - eoflgrn'o + lii - gimr, eo - (r) 

We display the corresponding contribution to the gi-F7r(gi) in Fig. 1 by the 
large-dashed curve. The solid line on this picture corresponds to the total 
result. 

Few comments are in order. As was expected, the interference term starts 
to grow at high enough momentum transfer. This is the main qualitative 
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effect. The contribution under consideration is subject to Sudakov correc- 
tions. An estimate of these corrections reveals that they are small enough. 
The reason for that is the large scale 45" which enter to the formula (0); thus 
the ^0 is far away from the end point region ±1 even at large q]_. As the 
second remark, we note, that the form factor is the integral characteristic 
and so, we expect that any reasonable regularization of 5 functions which 
respect constraints, will not change the F^^ drastically. 

At the same time we want to note that the result strongly depends on 
(^^). We display on Fig. 2 three different curves for q]_F-,^{q\^ which cor- 
respond to different choice of (^^). The small-dashed curve corresponds 
to the asymptotic wf with (^^) = 0.2 and g{ij) = 0. The large-dashed 
curve corresponds to (^^) = 0.3, g{fi) = 0.44. The solid line corresponds to 
(C^) = 0.4, g{fi) = 0.88. Let us note, that in these calculations we did not 
change the infrared regulator, which is an effective parameter for calcula- 
tion of the "standard" F^^ contribution. But we made all necessary substi- 
tutions related to changing of the parameter g{fi) in eq.(^). We observe the 
strong dependence on parameter (^^). In particular, at (^^) = 0.4, the soft 
contribution to the form factor at = lOGeV^ is equal to (l±FT^{q^) ~ 0.4. 
For the asymptotic wf with (^^) = 0.2 we have ^-F7r(^) — 0.1. These as- 
timations are in a very good agreement with absolutely independent recent 
calculations il^, who used a quite different method. 



4. Summary and Outlook. 

We hope we proposed the new "soft" mechanism which could be respon- 
sible for the explanation of the phenomenological success of the dimensional 
counting rules. The main idea is based on the very unusual feature (notice- 
able fluctuations of the transverse momentum ) of the wave function . It 
leads to the "two-hump" shape for the nonperturbative wf in the transverse 
direction and to the appearing of the new scale {S = O-SGeV"^) in the prob- 
lem, in addition to the standard low energy parameter {k\) ~ O.lGeV^ Q. As 
the consequence of it, the increasing of the in the pion form factor leads 
to the increasing of overlap of these humps. It gives the growing contribu- 
tion at intermediate ( parametrically it falls off very quickly, of course). 
Together with the standard, decreasing contribution, it could simulate the 
leading twist behavior in the extent region of g^. Our conclusions essentially 



support the picture described in [|18l, but from the quite different side. 

It would be interesting to check the conjecture formulated above (about 
an equivalence of the transverse momentum defined in terms of the gauge 

^ In the standard quantum mechanics it is very unhkely to have the two-particle lowest 
state wave function with such unusual form. However, we consider our wf as an effective 
one, which counts infinite number of soft gluons. The information about them is hidden 
into the definition of the effective transverse momentum. We believe that precisely these 
gluons are responsible for the making constituent quark from the current massless quarks 
and gluons. 
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invariant operatorors in QCD and perpendicular momentum in quark model) 
in the phenomenological analysis for the different amplitudes. Besides that, 
we would expect that the same feature might occur in the nucleon wf and 
we hope to discuss it somewhere else. I thank Tolya Radyushkin and Stan 
Brodsky for useful comments and insights on the number of subjects related 
to this letter. 

This work is supported by the Texas National Research Laboratory Com- 
mission under grant # RCFY 93-229. 
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FIGURES 



Fig.l Small-dashed line is the F^^ contribution to the pion form factor 
g^-F^^(g^). Large-dashed line is the interference term, ^F^^{q^). Solid 
hne is the total result ^F^{^) for gf(//) — 1. 

Fig. 2 Small-dashed line is the 9i-F7r(gi) at (^^) = 0.2; Large-dashed hne 
corresponds to (^^) = 0.3; Solid hne is the result for (^^) = 0.4. 
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